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This paper is a review of recent lattice determinations of the light quark masses. It describes the method 
employed to calculate quark masses in the lattice formulation, and the extrapolations required to reach the 
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1. INTRODUCTION 

The masses of the light quarks are fundamental 
parameters of the Standard Model, and yet are 
still surprisingly poorly determined. The Particle 
Data Book gives their values as 

m^(lGeU) = 2^8 MeV 

m^(lGeU) = 5 -> 15 MeV (1) 

TOl(lGeU) = 100 -> 300 MeV 

Over the last few years, several lattice projects 
have calculated these quark masses @-|. How- 
ever, these predictions have been mutually con- 
sistent only within a factor of two. It is therefore 
vital that the reason for this uncertainty is under- 
stood. This review describes the main reasons for 
the spread of values in terms of (i) the parame- 
ters used by the different collaborations, and, (ii) 
the extrapolations required to reach the physical 
regime from the region simulated. It is very rea- 
sonable to assume that the spread in lattice values 
can be greatly reduced in the next couple of years 
once more precise calculations are made, and the 
issue of the extrapolations is resolved. This would 
lead to extremely accurate lattice determinations 
of the light quark masses, far more constraining 
than eq.(Q) above. 

The bulk of lattice calculations have studied 
the isospin averaged quark mass mi = \(m u + 
rrid) and I will restrict my attention to these cal- 
culations. I note though that there has been 
a preliminary lattice investigation of the isospin 
mass splitting m c i — m u by using electromag- 
netic as well as gluonic fields. They obtain a value 
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for m u /md of « 0.5 in excellent agreement with 
other methods showing that isospin splitting ef- 
fects are also technically accessible on the lattice. 

In the next section I review the principles 
involved in lattice calculations of light quark 
masses. I then discuss the various extrapolations 
that are required in order to make contact with 
the real world in section 3. This section will be 
the main focus of the review. In section 4, I re- 
view the recent lattice results, attempting to ex- 
plain some of the inconsistencies that occur. In 
this section I give rough lattice values of the light 
quark masses, based on the data and the issues 
raised in section 3. In the last section I compare 
these values with estimates using other methods. 

For other reviews on this topic, including sum- 
maries of lattice results prior to 1994 see ||||. I 
note also that and [|llj contain very recent 
reviews of this topic. 

2. HOW THE LATTICE CALCULATES 
THE LIGHT QUARK MASSES 

The generic lattice approach to spectrum cal- 
culations studies two-point correlation functions, 
G2(t), of hadronic interpolating operators in a 
background sea of glue (and sea quarks in the 
case of "unquenched" calculations). These back- 
ground configurations are weighted by the appro- 
priate Boltzmann factor (which means that the 
approach is identical to calculation of correlation 
functions in statistical mechanics). Lattice calcu- 
lations are fully non-perturbative, and, in princi- 
ple, only involve approximations that can be sys- 
tematically improved. In practice the inclusion of 
dynamical sea quarks (i.e. virtual quark - anti- 
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quark pairs) is a technical headache because it in- 
volves a huge increase in computational require- 
ments. However, more and more lattice results 
are being published with dynamical sea quarks 
(also called "unquenching" ) and this trend is cer- 
tain to gain momentum. 

In the calculation of G%(t), if one were fortu- 
nate enough to have The Interpolating Operator 
of the hadron in question, then G2{t) would be 
The Two-point Function. However, when using a 
non-perfect operator, O, G 2 (t) receives contribu- 
tions from all hadronic states which have non-zero 
overlap with O. It is straightforward to show that 
C?2(i) has the following form 

G 2 (t)=Y J Z. l e~ M * 

i 

where the sum is over the hadronic states, and 
Mi and Zi are the hadronic mass and overlap 
of the i— th state. Note that since the calcula- 
tion is performed in Euclidean space-time, the 
excited states are exponentially suppressed with 
respect to the fundamental state (i.e. the expo- 
nentials have real arguments). This means that 
G%{t) becomes The Two-point Function only for 
sufficiently large values of t. Obviously the pa- 
rameters of the ground state, in particular the 
mass Mo = M, can be extracted by fitting (^(i) 
for t sufficiently large, i.e. 

G 2 {t)^Ze- Mt as t -> oo. 

Here comes the crucial point. This mass, M, is 
not the mass of the physical hadron. It is the 
mass of a hadron in a universe where the quark 
masses are the same as those input into the lattice 
calculation, i.e. the lattice is simulating a form of 
QCD, but not the QCD of the real world. 

In fact it is a little more complicated than this. 
The mass M (and all lattice predictions) are func- 
tions of all the input parameters of the lattice 
calculation. Specifically these parameters are: 

• quark masses, m; 

• lattice volume, V; 

• lattice spacing, a; 



• treatment of sea quark effects 
(i.e. quenched or unquenched); 

• type of lattice action used. 

In this section I consider only the dependence 
on the first quantity above; I leave the other ef- 
fects to the next section. 

Obviously we do not know the correct values 
of the quark masses to use in the calculation - 
these are, of course, the parameters that we wish 
to calculate. So the next step is to adjust these 
light quarks masses in the Lagrangian until the 
hadronic mass, M, agrees with the experimental 
values. Thus, for instance, we equate: 

My (m s ,mi) = M K , = 770MeV, (2) 

where My(m,, mi) is the lattice value of the vec- 
tor meson composed of quark with masses m s and 
mi . We are naturally free to use any hadron con- 
taining a strange quark in order to determine the 
strange quark mass, e.g.[] 

M v (m s ,m s ) = M 4> = 1020MeV, 
or 

M PS (m s ,mi) =M K = 494MeV, (3) 

where PS stands for pseudoscalar. The differ- 
ences in these various determinations of m s are 
systematic effects and should be included in the 
error determination. Note also that the deter- 
mination of m s from the K* and <f> are gener- 
ally in excellent agreement. This is due to the 
facts that (i) the experimental values of these 
hadronic masses are linear in their consitucnt 
quark masses, and (ii) to a very good approxima- 
tion, the lattice values of My (7711,7712) are also 
linear in mi and 7712. 

Normally, due to the computational expense, 
the lattice calculation is performed using only a 
small number (say 3 - 5) of quark masses. The ad- 
justment of the quark masses mentioned above is 
done by means of an interpolation/extrapolation 
of these data points. 

Note that there is an alternative approach to 
the calculation of light quark masses on the lat- 
tice using the axial Ward identity (see [Q & |l2| ) . 

2 Normally baryons are not used to set quark masses since 
their channels are more noisy in lattice calculations. 
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There has also been a very recent lattice recent 
calculation of the "renormalisation group invari- 
ant mass" in 

Once the values of the lattice quark masses 
have been found, they are then expressed in the 
MS scheme, traditionally at the scale 2 GeV. 
(For convenience, I note the conversion factor be- 
tween 2 GeV & 1 GeV is m(l GeV)/m(2 GeV) « 
1.4.) The matching procedure uses the relation- 
ship between the two schemes, 

m{p) = Z m (fj,a)m(a), (4) 

where m(p) is the quark mass in the MS scheme 
at the scale /i. Z m has been calculated in pertur- 
bation theory and its renormalisation group- 
improved form at next-to-leading order has been 
given in QJ^| Note that the expansion parameter 
g 2 relevant for the matching coefficient Z m is a 
short-distance (i.e. small) coupling constant, so 
perturbation theory should be fairly well behaved 
for this quantity (barring "tadpole" contributions 
|2| ) . This contrasts with the large distance scales 
governing hadronic quantities like M (which re- 
quire a non-perturbative tool such as the lattice 
to calculate correctly). 

There are two families of the lattice actions cur- 
rently used, Wilson and Staggered. While the 
staggered action has features that are particularly 
endearing for chiral physics, it has the big disad- 
vantage that its 0(g 2 ) coefficient in Z m is very 
large; typically this term is 50 - 100% of the zero 
order term ||. Therefore connecting any lattice 
quark mass value with continuum schemes intro- 
duces unknown higher order effects which are un- 
controlled. For this reason I do not consider re- 
sults from the staggered approach in this review. 

In the Wilson action, the quark mass is addi- 
tively renormalised. The connection between the 
quark mass parameter in the lattice Lagrangian, 
K , and the quark mass is 



3 Recently the lattice co mpo nent of Z m has been calcu- 
lated non-perturbatively [ha which means that the only 
perturbative component remaining in Z m is in the contin- 
uum matching. 



where K c corresponds to a massless pion, i.e. 
M 2 s (K c ,K c )=0. 

All fields and couplings in the lattice action 
have had their dimensions removed by multiplica- 
tion of appropriate powers of a. These powers of 
a must be put back in all dimensionful predictions 
(see eq.(|])). The lattice spacing a can be obtained 
by fixing the lattice prediction of any dimension- 
ful physical quantity with its experimental value. 
Thus one of the physical 'predictions' from the 
lattice is used to set a (this corresponds to set- 
ting the scale Kqcd), and one lattice 'prediction' 
is used for each quark mass determination (see 
above). The spacing a can be set from any di- 
mensionful physical quantity (that has been mea- 
sured by experiment!). Different physical quanti- 
ties can lead to slightly different values of a, and 
this systematic effect should again be reflected in 
the error estimate. 

In this section I have described, in principle, 
how lattice calculations of light quark masses are 
performed. In the next section, I review the de- 
tails of this approach, in particular the effect that 
the extrapolations in the quark masses, volume, 
lattice spacing etc. have on the final lattice pre- 
diction. 

3. EXTRAPOLATIONS 

3.1. Basic Ideas 

In the previous section I stressed that there are 
several input parameters in lattice calculations of 
strong inter ation physics: the quark masses, vol- 
ume, lattice spacing, type of lattice action used 
and the number of quark flavours included in the 
sea, Nf. 4 Assuming that all of these parameters 
were set to their physical value, then the lattice 
would give absolute predictions of QCD. How- 
ever, for technical reasons, these input parame- 
ters are almost never set equal to their physical 
values, and so the following extrapolations are re- 
quired (note that in lattice calculations, one is 
free to use sea quark masses, m sea , different to 
the "valence" quark masses, m v ): 

{m sea } -> {m u ,m d ,...} 

4 The quenched approximation corresponds to Nf = 0. 
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(G) 



Sea Quark Mass 



{m } -> {m u ,m d ,m s , ...} 

Volume — > oo 

a -> 

AT) -> "2". 

Thus the lattice simulates QCD, but with an 
unusual (i.e. unphysical) set of parameters - the 
light quark masses are larger than their physical 
values (typically they are i; 50 MeV) , the volume 
is finite, space-time is discretised, and the number 
of dynamical quark flavours does not always cor- 
respond with the real world. It is essential that 
all of the extrapolations eq.@ are under control 
in order for the lattice to make accurate determi- 
nations of physical quantities. 

A reliable extrapolation requires data points 
over a range of values of the parameter in ques- 
tion. Thus, for example, an extrapolation in a 
requires data at several values of a. Due to the 
present limitations in computing power, it is not 
possible to cover a large amount of the (multi- 
dimensional) parameter space in order to perform 
all of the extrapolations. Therefore most lattice 
projects concentrate on the effects of varying only 
one or two of the above parameters. Different 
groups often concentrate their efforts in studying 
and extrapolating different parameters and this 
fact can sometimes lead to different final predic- 
tions for the quark mass. The wide variation in 
lattice predictions for the quark masses referred 
to in the introduction is primarily due to this ef- 
fect. R Data at the same parameter values from 
different groups are generally compatible. 

Unfortunately, all of the extrapolations in 
eq.(||) are required in order to obtain The Lat- 
tice Prediction of the quark mass (or, in fact, of 
any physical quantity). Thus, while the lattice 
calculations performed so far are very encourag- 
ing, more work will be needed to reach this ulti- 
mate goal. Fortunately, many physical quantities 
aren't nearly as sensitive to these input parame- 
ters as the quark mass, and therefore don't suffer 
from these extrapolation effects to the same ex- 
tent. 

In the following subsections I will discuss the 

5 In my opinion this point has not been fully communicated 
to the non-lattice community. 
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Figure 1. Schematic plot of the volume effect 
when naively extrapolating in sea quark mass in 
unquenched simulations. 



effects of extrapolations in the lattice spacing, a, 
the sea quark mass, m sea , as well as the related 
systematics due to unquenching effects. I do not 
consider volume effects, since they are relatively 
small. The extrapolations in m v are basically 
those discussed in the previous section. 

3.2. Lattice Spacing Effects 

In ||, a collection of world data for rni(2GeV) 
and rn^(2GeV) was extrapolated to the a = 
limit. They found that the quenched Wilson data 
had a strong dependence on a. Assuming a linear 
dependence m(a) ~ 757(0) x (1 + xa), the contin- 
uum extrapolant, 777(0), was as much as 40% lower 
than the simulated data. (This applied to both 
777/ and 777 s .) An analysis of unquenched Wil- 
son and staggered data was presented as "prelim- 
inary" since these simulations were in their early 
stages j§. 

In J|]7 some evidence for an 0(a) effect was 
found and this was included their final quark mass 
estimate. 

Another group, ||, analysed quenched Wilson 
and improved Wilson m s data from the APE col- 
laboration over a range of fairly fine lattice spac- 
ing. They found that any 0(a) effects were at 
least as small as the statistical errors in the cal- 
culation and therefore no extrapolation in a was 
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performed. 

It is my belief that more data of better quality 
is required to unambiguously determine the quark 
mass dependency on a. 

3.3. Sea Quark and Unquenching Effects 

A very interesting study of the dependency of 
mi and m s on the sea quark mass was presented 
in H . The calculation was at a fixed value of the 
strong coupling constant (naively corresponding 
to a fixed value of a) and used three values of 
the sea quark mass, the heaviest being at around 
the strange quark mass, and the lightest being 
approximately half m s . A parallel calculation was 
performed at a similar value of a in the quenched 
approximation (i.e. corresponding to Nf = 0, 
or equivalently, m sea — oo). The dependence of 
physical quantities such as meson masses on m sea 
was found to be strikingly large. This implies that 
the lattice spacing obtained, for instance, from 
M p (see Sect. ||) is itself a function of m seo .^ 

A little thought shows that this dependence of 
a on m sea should not be surprising. It has long 
been known that, in order to maintain roughly 
fixed physics (i.e. constant a values) when com- 
paring quenched and unquenched data, the cou- 
pling (3 — G/g 2 had to be renormalised by, say, 
~ 0.6 (see e.g. Jill ). This implies that, on the 
m sea — (3 plane, there is a curve of constant a 
connecting the quenched (m sea = oo) simula- 
tion at, say, /3 = 6.0 with the unquenched point 
(corresponding to a small value of m sea ) with 
(3 « 6.0 — 0.6. Now consider varying m sea while 
keeping [3 fixed. This is the parameter space in- 
vestigated by the SES AM collaboration || . This 
sweeps along a line with a component orthogonal 
to the curve of constant a. Therefore it is en- 
tirely reasonable that one finds a variation of a 
along this direction as reported by || . 

This feature of a variation of a with m sea has 
interesting implications. It means that as one ex- 
trapolates in m sea to the physical value, mi, the 
lattice volume decreases in physical units (since a 
decreases) . This is displayed schematically as the 
dashed curve in Fig.l. Normally one uses PCAC 



6 The strong dependence on m sea on the static quark po- 
tential and hadron masses has als o recently been observed 
by the UKQCD collaboration MM. 



or quark model considerations to model the be- 
haviour of masses with mi - at fixed volume which 
would correspond to solid vertical line in this fig- 
ure. In principle, this procedure will have to be 
modified to compensate for the finite volume ef- 
fects which mix with the m sea effects. However, 
in this relatively early stage of unquenched cal- 
culations, one is forced to ignore such effects in 
order to make progress. 

The SESAM collaboration also has data at non- 
degenerate "valence" quark masses m\ , m\ M. 
This means that they are able to apply eqs.(g& 
[}]) directly to obtain m s rather than relying on 
SU(3) flavour symmetry (see, e.g. p7|). Impor- 
tantly also they have data for m v ^ m sea and 
they present arguments about the correct proce- 
dure for extrapolating in m v and m sea pointing 
out the perils of following an incorrect procedure 
||. I reproduce their arguments here. 

Their explanation can be understood by study- 
ing Fig. 2 which is taken from Fig. 4 in the firsr 
reference in || (and slightly generalised). In this 
plot, the sea and valence lattice quark mass pa- 
rameters 1/K V & 1/K sea (see eq.(^)) are plot- 
ted against each other. The solid diagonal line 
corresponding to K = K sea corresponds to the 
degenerate case where much of the earlier sim- 
ulations have been performed. The three diag- 
onal dashed lines correspond to the constraints 
M PS (K v ,K sea ) = 0, M v (K v ,K sea 
M v (K v ,K sea ) = Mt/, (labelled |T|, |T| & [I 
spectively). Note that, as constructed in Q, these 
lines are parallel, and empirically their slopes are 
« —1. The points labelled u p" and "</>" are the 
physical points; p consisting of two light quarks in 
a sea of light quarks, and the 4> consisting of two 
strange quarks in a sea of light quarks. The point 
marked "0" corresponds to m v — m sea = 0. 

The two displacements on the y-axis denote the 
correctly determined quark masses mi and m s (see 
eq.(||)) - i.e. they are calculated in a sea con- 
sisting of physical "Z" quarks. The vertical dis- 
placements in the body of the figure show the 
incorrectly determined values of mi and m s . The 
incorrect mi value, marked "mP", is obtained if 
one performed a "quenched-like" analysis - i.e. an 
extrapolation in m v at fixed m sea . Note that this 
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Figure 2. A slightly generalised version of Fig. 4 
in [6] in which the sea and valence lattice quark 
mass parameters 1/K V & 1/K sea (see eq.(||)) are 
plotted. The three diagonal dashed lines corre- 
spond to the constraints M PS (K V , K sea ) = 
(denoted |T|) 
M v (K v ,K sea 



) = M p , ( 2 ), and 



( 3 ). The physical points 



M v (K l 

are labelled p and <f> with the corresponding cor- 
rect values of mi and m s . "O" corresponds to 
m v = m sea = ^ gee Sect.O for details. 



is around twice the correct value of m/. The mj 
value is the corresponding "quenched-like" value 
of m s . Note that this is surprisingly close to the 
correct m s value ! Finally, m d s ea is the m s value 
obtained if one worked with degenerate sea and 
valence quarks. This can be seen to be roughly 
half of the correct m s value. 

3.4. Summary of Extrapolation Effects 

In concluding this section, it is clear that fur- 
ther work is required to fully understand and 
compensate for the effects of the extrapolations 
of the input parameters to their physical values 
in a lattice simulation. As detailed above, both 
the continuum extrapolation (a — » 0, see Sect. 
IjO) and the sea quark effects (Sect. 3.3) intro- 



that are the main source of the large spread of lat- 
tice values for the light quark masses presented in 
contemporary literature. 

Despite these uncertainty, lattice calculations 
performed by different groups at common param- 
eter values are in generally excellent agreement 
leading to great optimism for the future. 

4. RECENT LATTICE RESULTS 

For the reasons outlined in Sect. ||, I will re- 
strict myself to discussing recent lattice results 
from the Wilson family of actions. There have 
been many calculations of light quark masses over 
the years (see |^,^| and references therein), but I 
will concentrate on the more recent results. These 
come from 4 groups [|-| and are shown in Ta- 
ble [|. (Note that the APE/Rome group has new 
data in |l2| which are compatible with those from 
&|[.) Alongside the values for Wl(2GeV) and 
m s (2GeV) are indicated whether (a) the simula- 
tions were performed in the quenched approxima- 
tion, (b) the correct m sea extrapolation was per- 
formed (see Sect. 3. 3) (c) an attempt was made 
to allow for 0(a) errors (see Sect.[3.2D, and, (d) 
whether an "improved action" was used[|. Ide- 
ally, the answers to these questions should all be 
"Yes" , but, as can be seen, there is no one simu- 
lation w here this is the case. (See the discussion 
in Sect. 



The quotation marks around the unquenched 
results for the LANL and FNAL groups signify 
that these are considered as preliminary |ft[, or 
are based on early unquenched results and plau- 
sibility arguments [Qj. 

It is interesting to note that the small value for 
the unquenched mj(2GeV r ) in || and Q arc not 
reproduced in the fuller unquenched analysis of 
||. This is due to two extrapolation systemat- 
ics: (i) the sea quark effects are correctly taken 
account of in , but not in || (see the discus- 
sion surrounding "mf e9 " in Sect. 3.3); and (ii) 
the 0(a) effects are treated differently - Q does 
no a — > extrapolation since only one value of 
the lattice spacing was used (furthermore an tin- 
improved action was used), whereas both p] and 



duce potentially large systematics that can be as 
much as a factor of two. It is these extrapolations 



7 These are lattice actions with reduced discretisation 
systematics. 
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Table 1 

Recent lattice results for the light quark masses. 



Group 



mj(2GeV) 
[MeV] 



rrT s (2GeV) 
[MeV] 



Unquenchcd 



Correct m s 







extrapolation ? taken ? 



•'Improved Action" 
used ? 



APE [2] 
APE [5| 
LANL || 

FNAL g 

SESAMil 



3.4(4)(3) 

"2.7(3)(3)" 

3.6(6) 

"2.1^3.5" 

5.5(5) 

2.7(2) 



128(18) 

122(20) 

100(21)(10) 

"68(12)(7)" 

95(16) 

"54^92" 

166(15) 

140(20) 



No 
No 
No 
Yes 
No 
Yes 
No 
Yes 



No 
No 
Yes 



"Yes" 
Yes 
Yes 

"Yes" 

No 
No 



Yes 
Yes 
Yes& No 

Yes 

No 
No 



Table 2 

Recent non-lattice results for mi and m s 
QCD Sum Rules 



[18] 
[19] 
J1S[ 
[20| 

[21] 
[22] 
J23] 



rh~i(2GeV) = 4.9(l.)MeV 
3AMeV < m(2GeV) < 5.6MeV 
m(2GeV) 5MeV 



m;{2GeV) = 114(21)MeV 
w;{2GeV) = 141(21)MeV 
r?r;{2GeV) = 146 (14) MeV 
m;{2GeV) = 102(13)MeF 



Chiral Perturbation Theory 



[24] ms/mi = 24.4(1.5)Mey 



Experiment 



[25] m;(2GeV) « imMeV 



|4| incorporate a linear a — > extrapolation pro- 
cedure (see Sect. 3.2). 

The other interesting feature to note is that 
the relative quenching effects in wn:(2GeV) ap- 
pear small, whereas they appear significant in the 
case of mi(2GeV) || . This can again be under- 
stood in terms of the discussion in Sect. 3.2 sur- 
rounding Fig. 2. 

By considering the extrapolation effects dis- 
cussed in Sect. 3, I give the following rough esti- 
mates for the values of the (unquenched !) quark 
masses based on the results in Table |l|. 



m^GeV) = 2^5 GeV, 
mT s (2GeV) = 80 -> 140 GeV. 



(7) 
(8) 



the values in Table [l]. The rn^(2GeV) estimate is 
obtained by considering the fact that the LANL 
(unquenched) predictions are likely to be too low 
due to their sea quark treatment, whereas the 
SESAM value is possibly too high due to 0(a) 
effects. The values in eqs.(0,||) are obviously only 
meant as a guide, and until the parameter space 
of lattice calculations is more fully explored, it 
is not possible to perform all the extrapolations 
required to determine The Lattice Prediction of 
mi{2GeV) and m;{2GeV). 



5. COMPARISONS 
METHODS 



WITH 



OTHER 



The mi(2GeV) range is chosed so that it covers 



I give in Table ^ some determinations of light 
quark masses from non-lattice methods. (These 
have all been run to the scale 2 GeV.) Note the 
preliminary result from the Aleph Collaboration 
for m s , |25). The lattice estimates in eqs.(^fe^) 
are in fairly good agreement with these results, 
although they seem generally on the low side. 
This point has been investigated in |^6| where 
lower bounds on mi and m s are derived from 
fairly fundamental principles. The lower limit for 
rh~i(2GeV) in eq.(0) is in some conflict with the 
bounds in @. (See also |7|.) 
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DISCUSSIONS 

A. P. Contogouris, University of Athens 
Some time ago, P. Lepage of Cornell Univ. 
proposed a method for lattice calculations and 
claimed that it could be used in a small computer, 
even a laptop(l). Did you use his method? And 
if not, why? More generally, would you like to 
comment about the present stage of Lepage 's ap- 
proach? 
C.R. Allton 

No, the calculations that I have described do not 
use Lepage's "coarse-lattice" method. This is a 
proposal that adds extra terms to the lattice ac- 
tion in order to reduce the effects of discretisation 
- i.e. it is a so-called "improvement" scheme. 
The important feature of Lepage's approach is 
that the coefficients of these additional terms are 
set by a re- arranged perturbation series (contrast- 
ing with other improvement methods which use 
a fully non-perturbative technique). Lepage's ap- 
proach has had success in reducing the rotational 
non-invariance of the static quark potential, and 
in quantities heavily dependent on that potential 
(such as physics of the heavy sector). However, 
progress in the light quark sector has been limited 
so far and this is the reason why the method has 
not been used for the determination of light quark 
masses. 

K.F. Liu, University of Kentucky 
/ want to add to what Chris has said. While 
the improved action program has been working re- 
markably well for the pure gauge sector and heavy 
quarks for NRQCD, its progress for light quarks 
is still limited at this stage. 



